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Abstract. We would like to formulate relativistic dissipative hydrodynamics for multi- 
component systems with multiple conserved currents. This is important for analyses of the hot 
matter created in relativistic heavy ion collisions because particle creations and annihilations of 
various particle species are frequently taking place there. We show that consistent formulation 
in such systems involves many non-trivialities, and derive constitutive equations that satisfy 
Onsager reciprocal relations and describe the systems without ambiguity. 



1. Introduction 

One of the most important discoveries of Relativistic Heavy Ion Collider experiments is that 
relativistic ideal hydrodynamic models quantitatively describe the collective behavior of the 
quark-gluon plasma (QGP) created in the collisions [1J. The success is a great motivation to 
take non-equilibrium processes into account to explain small discrepancies between theoretical 
predictions and experimental data. Large Hadron Collider experiments also expect heavy ion 
collisions of higher energies, and viscous hydrodynamic models would become more important 
as one might observe less strongly-coupled QGP because of asymptotic freedom in QCD. 

The ambiguities of viscous hydrodynamic models lie in the formalism of relativistic dissipative 
hydrodynamics itself. If one considers only the first order dissipative correction to the entropy 
current, theories are known to become acausal and unstable [2]. Thus varieties of theories with 
the second order correction have been proposed. However, the resulting constitutive equations 
for the dissipative currents differ among theories. Also, most viscous hydrodynamic models are 
based on the theories for single component systems with binary collisions. Considering that 
the QGP is a multi-component system with particle creation and annihilation, it is important 
that one establishes a multi-component formalism which takes account of inelastic scatterings. 
This naturally motivates us to consider multiple conserved currents because one needs to solve 
continuity equations for conserved charges, rather than the ones for particle numbers. 

2. Formulation of Relativistic Dissipative Hydrodynamics 

We would like to formulate relativistic dissipative hydrodynamics for multi-component and 
multi-conserved current systems [3\- Relativistic hydrodynamics is a general theoretical 
framework which describes collective motion of strongly-coupled systems with the energy- 
momentum conservation d^T^ v = 0, the charge conservations d^Nj = and the law of increasing 



entropy d^s^ > 0. Here the index J denotes different types of conserved quantum number. 
Typically, one can consider the baryon number B or the strangeness S for such quantities in 
relativistic heavy ion collisions. 

We introduce tensor decompositions of the energy-momentum tensor T^ u and the conserved 
currents Nj in terms of the flow field to define thermodynamic quantities 

T» v = (e + Se)uV - (P + n)A"" + W »u v + W v u» + tt^, (1) 
N$ = (n J0 + 5nj)u» + V?, (2) 

where A^ = g^ u — u^u v is the projection operator. There are 2 + N equilibrium quantities: 
energy density eo, hydrostatic pressure Pq and type-J charge density tijq. On the other hand, 
one has 10 + AN dissipative currents: bulk pressure II, energy density distortion Se, energy 
current W^, shear stress tensor ■K tJ,v , type-J charge density distortion Snj and type-J charge 
current Vj. Although thermodynamic stability conditions [HE] require Se and Snj to be zero, 
we have to keep them for the moment to correctly count the number of the dissipative currents. 
In ideal hydrodynamics, the conservation laws and the equation of state Po = Po(eoi { n Jo}) give 
sufficient number of equations to describe systems because the number of unknowns is 5 + N in 
such systems. In dissipative hydrodynamics, one needs constitutive equations because there are 
the 10 + 4iV additional unknowns II, Se, W^, tt^ u , 5nj and Vj. We would like to derive them 
from the law of increasing entropy, because it embodies irreversible processes. 

The energy-momentum tensor and the conserved currents are expressed in relativistic kinetic 
theory with the phase space distribution /* as 

where g\ and qf are degeneracy and type-J quantum number of i-th particle species, respectively. 
We introduce the higher moment equations for the energy-momentum conservation and for the 
charge conservations 

da I»™ = Y,J -0^Mp?d a f=Y^, (4) 

= Y.j^im^f=y^ (5) 

to derive second order constitutive equations. Y^ v and Yj are to be determined from the law of 
increasing entropy. The conventional Israel-Stewart theory [5] has only Eq. @ because single 
component systems with binary collisions are considered. However, one has to introduce the 
new sets of equations ([SJ to consistently match the numbers of equations with that of dissipative 
currents in multi-component and/or multiple conserved current systems. 
The entropy production can be written in kinetic theory as 



V = -a,E / ^Ir^f/ 4 m/* -^(1 + ^)111(1 + ^)] =£/ ^wfiM (6) 



(2tt)3^ 



(2vr)3^ J 



where y l is defined in /' = [exp(y l ) — e] _1 . The sign factor e is +1 for fermions and —1 for 
bosons. We estimate the distortion of distribution 5p = f % — /q in terms of dissipative currents 
by generalizing the Grad's moment method [S] for the systems with multiple conserved currents. 
If we require that the resulting constitutive equations satisfy Onsager reciprocal relations [6l[7], 
the only possible moment expansion would be 

Sf = -/5(l + e/5)(pj , E^+^ e /»'). ( ? ) 

J 



where and e^ v are self-consistently determined so that the distribution reproduces the 
energy- momentum tensor and the conserved currents ([3]). Then the e's are expressed as linear 
combinations of the dissipative currents. On the other hand, the entropy production Q is 
expressed in terms of e's and Y's as 

V = T, £ tYj+ev P Y vp >0. (8) 
J 

The law of increasing entropy requires that y's are linear combinations of e's. Then the second 
order constitutive equations for multi-component systems with multiple conserved currents are 
obtained by writing the moment equations (HJ) and ([5]) with the dissipative current: 

n = -(V^-CnseD^+YsCnsnjD^ 

j 

- mDU + £ x&ILD^ + Xim^ + xfiu^^ 
J 

+ E + xW^V^ + xW^«" + XrWV^ 

+ E x&jVj^ff + e xVX v 4 + e T&vjVjlDvr + E 4n d 

J,A J J J 

+ Xn^vK^, (9) 



J 



twA^DW v + E + X b wwW»D± 



+ XW^V/ + + XwwW^u, 

- £ iWvA-'DVj + E xtffc + E xWr^ 

J J,K J 

+ E x c WVj vfvuv + E xW/Jvy + E XW^J 

J J J 

J 

+ ^^nv^ + xWnv^ + xWn^ + xwnV^n, (10) 



e t VjVk a^dv u k + e xtv K v%D^ + x: 

A' A,L A' 



E xWY£V"U" + E xW V^"" + E ^ v^Vt 

AAA 



rp ' /vyjvv rp 

A 



+ Xv jW W^ u u u + Xv jW W v V v u» + Xy jW m 

+ E X&^V^ + X^^V^ + X^*^ + X^A^V^TT.p 
A 



E xghnv"^ + x^nllV^ + x ^ n iW + xfcnV^H, (11) 



if 



J 

+ E X a / Vj Vj^ ff + E X^V"> i + E X^>^ + E xtv^V? 

J,K J J J 

+ x^rnnv^u"). (12) 



Here D = ii^<9^ and = A^ v d v are time- and space-like derivatives. The angle brackets denote 
the operation A^Bfi = [^A^A^ + A^A^J - iA^A^A"^. C's, «'s and rj are first order 
transport coefficients; C is bulk viscosity, kww energy conductivity, KVjVj charge conductivity 
and r] shear viscosity, r's are relaxation times, and x' s are second order transport coefficients. 
The ratios of the first and the second order coefficients can be determined in kinetic theory. 

3. Discussion and Conclusions 

We have derived constitutive equations for multi-component systems with multiple conserved 
currents in arbitrary frame. We would like to emphasize four points. Firstly, the equations 
have Onsager cross terms, although most of the conventional formalisms do not consider them 
carefully. These terms contain important physics such as Soret and Dufour effects, which are 
chemical diffusion caused by thermal gradient and its counter-effect respectively. The reciprocal 
relations are satisfied. Secondly, there are relaxation terms - terms with time-like derivatives of 
the dissipative currents themselves - in the constitutive equations. These terms play important 
roles in preserving causality because they cause exponential damping of the dissipative currents. 
Thirdly, we have independent equations for energy and charge currents. This means that our 
formalism can describe systems in an arbitrary frame unlike the conventional formalisms in 
which one has to choose either energy or particle frames. Fourthly, we have several second order 
terms which do not appear in Ref . [5] . For further detailed discussion including comparison with 
other approaches [H [91 [10] , see Ref. [5] . 
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